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基于 双 伸 缩 因子 映射 的 小 波 分 形 图 像 压 缩编 码 方法 


惠 存 阳 
(西安 文理 学 院 数学 系 ， 西 安 710065) 
摘 要 : 本 文 基于 小 波 变换 研究 了 分 形 图 像 压缩 系统 .根据 小 波 系数 的 能 量 分 布 特性 ， 给 出 了 新 的 小 波 树 
定义 与 分 类 方法 ， 并 在 小 波 域内 建立 了 具有 双 伸 缩 因子 的 压缩 映射 . 在 将 新 的 压缩 映射 应 用 于 图 
像 编码 过 程 时 ， 结 合 小 波 嵌 入 式 零 树 编码 思想 给 出 了 新 的 小 波 分 形 图 像 压 缩 方法 .数值 实验 给 出 


了 本 文 方法 与 已 有 方法 的 比较 . 
关键 词 : 分 形 图 像 编码 ， 伸 缩 因子 ， 小 波 ;， 压缩 映射 
分 类 号 : AMS(2000) 42C40 中 图 分 类 号 : 0184; 0174.3 文献 标识 码 : A 
1 引言 


分 形 图 像 编 码 是 基于 压缩 映射 的 一 种 图 像 编 码 方法 ， 由 Barnsleyi 首先 提出 .该 方法 常常 
能 够 达到 较 高 的 压缩 比 ， 优 于 其 它 编码 方法 . 但 Barnsley 并 未 公布 其 具体 算法 . 之 后 ， 他 的 学 
^E Jacquin?) 基于 图 像 分 块 的 方式 ， 采 用 逐 块 变换 系统 (PTS)， 首 次 给 出 了 分 形 图 像 编码 的 算 
ik. 但 应 用 PTS 时 ， 需 要 进行 大 量 的 搜索 计算 ， 耗 时 较 多 . 为 此 ， 有 许多 文献 给 予 了 改进 ， 
如 为 了 减 小 搜索 计算 量 而 提出 的 固定 网 络 搜索 法 、 分 类 搜索 法 、 四 又 树 搜索 法 以 及 最 优 搜索 法 
等 。 在 子 块 划分 方面 ，Fisher ADS] 提出 了 自 适应 四 又 树 编码 方法 ， 这 就 是 目前 常用 的 分 形 图 
形 压缩 方法 . 进一步 的 工作 可 参见 文献 [4-8]. 

分 形 图 像 压 缩 方 法 发 展 的 另 一 个 重要 方向 是 结合 小 波 分 析 的 编码 方法 研究 ， Davislg] 首先 讨 
论 了 小 波 域 内 的 分 形 图 像 编 码 方法 ， 并 对 分 形 编码 方法 和 解码 收敛 性 进行 了 研究 ， 该 方法 较 空 
闻 域 的 分 形 编 码 方法 有 了 许多 改进 .首先 是 编码 时 间 缩 短 ， 其 次 通过 适当 的 小 波 选择 可 使 恢复 
图 像 质 量 没有 方块 效应 ， 主 观 上 图 像 具 有 较 高 的 质量 .但 Davis 所 采用 的 压缩 映射 依赖 于 单 伸 
缩 因子 ， 给 匹配 造成 困难 ， 常 常 需要 分 解 子 树 ， 使 得 计算 量 增加 . 

我 们 注意 到 ， 自 然 界 中 的 图 像 并 不 是 数学 意义 上 的 严格 自 相 似 图 像 ， 因 而 具有 单 伸缩 因子 
的 仿 射 变换 在 实际 图 像 压 缩 中 常常 会 产生 较 大 的 匹配 误差 .为 了 提高 子 树 与 父 树 的 匹配 精度 ， 
减 小 巨 配 搜索 范围 ， 本 文 根 据 小 波 性 质 ， 研 究 了 新 的 小 波 树 分 类 方法 与 双 因 子 匹 配 映 射 ， 实验 
表明 ， 本 文 所 给 出 的 双 因 子 匹配 映射 能 够 明显 的 减少 子 树 与 父 树 的 匹配 误差 ， 在 给 定 匹配 误差 
时 有 更 多 的 子 树 能 找到 匹配 父 树 ， 因 而 有 利于 提高 恢复 图 像 的 质量 及 编码 速度 . 





2 Davis 小波 分 形 图 像 压缩 方法 


2.1 小 波 树 

Davis 将 三 个 方向 上 的 小 波 系数 定义 为 一 棵 小 波 树 ， 如 图 1 所 示 . 在 四 级 小 波 分 解 时 ， 小 
波 树 有 三 个 根 结 点 均 位 于 第 三 层 ， 每 个 根 结 点 在 较 细 尺 度 的 同一 方向 上 又 有 四 个 “孩子 ” ， 
图 1 中 用 方 格 状 块 表示 所 有 的 根 结 点 与 父 结 点 . 

收 稿 日 期 : 2009-05-18. 作者 简介 : BARWNE) B, WAR. 研究 方向 : 小 波 分 析 及 其 应 用 . 
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与 空间 域 的 子 块 及 父 块 概念 相对 应 ， 在 小 波 域 有 子 树 (range subtree) 与 父 树 (domain sub- 
tree) 的 概念 .与 子 树 比较 ， 父 树 的 根 结 点 位 于 尺度 较 粗 的 分 解 层 上 . 如 图 1 所 示 ， 若 根 结 点 
EJ = 3 的 分 解 层 上 的 小 波 树 称 为 子 树 ， 而 根 结 点 在 分 解 层 了 = 4 上 的 小 波 树 称 为 父 树 . 所 有 
父 树 构成 的 集合 称 为 父 树 库 . 








图 1: 小波 树 的 构造 


2.2 ”压缩 映射 与 图 像 压 缩 思想 

由 于 小 波 树 的 特殊 结构 ， 对 父 树 的 抽样 及 仿 射 变换 与 空间 域 有 所 不 同 ， 在 小 波 域 中 ，““ 平 
均一 抽样 ”操作 4 采用 减少 一 个 分 辨 率 的 办 法 进行 ， 把 父 树 的 “叶子 ”， 即 最 精细 尺度 上 的 小 
波 系数 置 为 0， 其 余 尺 度 上 的 小 波 系数 移 向 下 一 个 精细 尺度 . 经 过 4A 操作 后 ， 父 树 的 大 小 与 子 
树 相 同 . 

小 波 域内 的 仿 射 变换 与 空间 域 中 的 8 种 仿 射 变换 四 直接 对 应 .因而 旋转 及 反射 变换 涉及 到 
小 波 树 的 三 个 方向 ， 除 了 对 每 一 个 方向 上 小 波 树 的 子 块 进行 反射 、 旋 转变 换 外 ， 对 某 些 变换 还 
需要 交换 LH 与 HL 方向 上 子 块 的 位 置 ， 这 些 变换 称 为 等 距 算 子 ， 记 为 L， 小 波 域 中 的 质量 映 
射 定义 为 : (4D) = aAD， 其 中 a 称 为 伸缩 因子 ，DD 为 父 树 . 综 上 所 述 ， 小 波 父 树 到 子 树 的 
仿 射 变换 可 表示 为 

t(D) = QoLoA(D)= a(A(D)), (1) 
小 波 分 形 图 像 压 缩 的 基本 思想 是 ， 对 于 给 定 的 子 树 尺 ， 用 父 树 来 逼近 子 树 ， 即 寻找 适当 的 伸缩 
因子 、 等 距 变换 地， 使 
lg — aL(Â(D))|| 


达到 最 小 . 


3 ”小波 分 形 图 像 编码 方法 的 改进 


3.1 小 波 树 的 分 类 

我 们 知道 ， 在 对 图 像 作 小 波 塔 型 分 解 后 ， 三 个 方向 上 的 能 量 分 布 是 不 均衡 的 ，HH 方向 分 布 
的 能 量 最 小 ， 而 LH 方向 与 HL 方向 分 布 的 能 量 较 大 . 为 此 ， 我 们 定义 每 个 方向 上 的 小 波 系数 
为 一 棵 小 波 树 ， 并 将 小 波 树 分 为 如 下 四 类 : 
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1) 零 数 ， 相 对 于 阐 值 TT 树 上 每 个 系数 都 是 无 效 的 ; 
2) ” 根 重 要 树 。 只 有 根 结 点 是 有 效 的 ， 而 其 他 结 点 是 无 效 的 ; 
3) 可 自 变 换 树 : 如 果 它 可 找到 满足 匹配 误差 要 求 的 父 树 ; 

4) ”复杂 树 : 除 上 述 三 类 树 之 外 的 小 波 树 . 

5 EZWUO 的 分 类 方法 相 比 ， 本 文 分 类 方法 引入 了 根 重要 树 . 

在 编码 时 ， 最 低频 的 逼近 系数 量化 后 直接 保存 .小波 树 扫描 从 最 低频 ( 子 带 ) 开始 按 之 字 型 
进行 ， 每 一 子 带 内 的 位 置 扫描 完 后 进入 下 一 子 带 扫描 ， 直 到 最 高 频 CT) 为 止 ， 而 在 每 一 分 解 
层 内 ， 先 扫描 HL 方向 的 根 结 点 ， 再 扫描 世英 方向 的 根 结 点 ， 最 后 扫描 HH 方向 的 根 结 点 . 

3.2 ” 父 树 到 子 树 的 映射 

我 们 定义 父 树 到 子 树 的 映射 为 如 下 形式 的 仿 射 变 换 


T = diag(ki,ko,--- kn) o L, (2) 


Jrb n [CK BUR T-T- BUS. MET VERSOS I eM k 需要 的 比特 数 ， 本 文 考虑 两 个 伸 绾 因 
子 的 情况 ， 即 前 1 个 取 同 一 个 值 ， 后 -1 个 取 同 一 个 值 . 

记 父 树 为 D = (D'D?), THAR = (R, RP), PR RKA ACD?) 的 大 小 相同 ， 则 父 
树 与 子 树 的 匹配 误差 为 


1 n 
MSE = y ET tegl ks, ski; 2， K)LAD[ 


x; (IRI? — 21 G8, LÀ(D)) + 1D P — 2ka(R?, L(Â(D?))) + KLAQX?), (3) 
EFI ACGREXIGBES, M = dim(R). 

3.3 ”编码 与 解码 过 程 

根据 上 述 小 波 树 分 类 方法 ， 我 们 结合 嵌入 式 零 树 编码 加 思想 给 出 了 如 下 分 形 编码 方法 ， 假 
设 判别 小 波 系数 是 否 有 效 的 阔 值 为 九 ， 匹 配 误 差 闪 值 ， 则 有 具体 算法 如 下 : 

步骤 1 对 原 图 像 作 N 级 小 波 分 解 ; 

步骤 2 ”量化 并 保留 逼近 系数 与 最 粗 尺 度 上 的 小 波 系数 ， 

步骤 3 分 别 构造 子 树 库 与 父 树 库 ， 置 dep = N 1; 

步骤 4 扫描 第 dep 层 的 小 波 系数 ， 记 根 结 点 为 2 的 小 波 子 树 为 RS. RE CARRIE YN 
后 辈 ， 则 不 编码 ， 按 顺序 考虑 下 一 子 树 ， 否 则 对 R, 做 如 下 处 理 : 

1) 如果, 是 零 树 ， 则 记 下 其 类 别 ， 作 为 编码 ， 按 顺序 扫描 下 一 子 树 ; 

2) 如 果 ,是 根 重要 树 ， 则 记 下 其 类 别 ， 并 保留 根 结 点 ， 编 码 完成 后 ， 再 按 顺 序 扫描 下 一 
子 树 ; 

3) 如果 忆 ,是 可 自 变 换 的 子 树 ， 则 保留 父 树 D, 的 位 置 ， 两 个 伸缩 因子 及 仿 射 变换 类 型 作 
为 子 树 R 的 编码 ， 青 按 顺 序 考虑 下 一 子 树 ; 

4) 如 果子 树 是 复杂 树 ， 则 保留 根 结 点 ， 并 将 子 树 Rp 留 到 下 一 个 较 高 频 子 带 进行 编 
码 . 重复 1) 到 4) 操作 ， 直 至 第 dep 层 内 的 子 树 处 理 完 为 止 ; 

步骤 5 如 果 dep = 1， 编 码 结束 .和 否则， 置 dep = dep -- 1， 转 到 步骤 4. 

小 波 分 形 图 像 编码 流程 图 如 图 2 所 示 . 

解码 过 程 如 下 : 

步骤 1 任 取 一 图 像 TU， 将 保留 的 逼近 系数 与 小 波 系数 置 于 左上 角 ; 

步骤 2 重复 如 下 操作 : 
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1) 造 父 树 库 ; 2) HER TU 的 相应 子 树 按 其 分 类 进行 解码 ; 
直至 满足 要 求 为 止 ; 
步骤 3 ”做 小 波 逆 变换 ， 给 出 解码 图 像 . 








作 N 级 小 波 分 解 , 置 dep=N-1 


按 之 字 型 扫描 dep 层 小 波 系数 ,构造 小 波 树 


图 2: 改进 后 的 小 波 分 形 图 像 编码 流程 图 


4 实验 结果 


小 波 基 函 数 选 为 B9/7 双 正 交 小 波 Li1， 原 始 图 像 为 512 x 512Lena 头 像 ， 灰 度 值 量化 
为 256 级 ， 如 图 3 所 示 . 在 以 下 各 实验 中 均 作 6 级 小 波 分 解 . 

4.1 编码 方法 的 效果 

在 实验 中 ， 当 子 树 根 结 点 位 于 J = 2 的 分 辩 率 层 时 (参看 图 1)， 不 再 用 父 树 编 码 子 树 ， 而 是 
直接 编码 系数 . 父 树 位 置 与 仿 射 变换 类 型 用 等 长 编码 ， 如 对 第 五 层 的 子 树 ， 编 码 其 父 树 的 位 置 
需要 6bit， 下 一 层 需要 8bit， 储 存 每 一 仿 射 变换 用 3bit. 对 于 子 树 的 类 型 及 保留 的 小 波 系数 ( 先 
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量化 ) 用 算术 编码 进行 编码 . 实验 结果 如 表 1， 图 4 所 示 ， 其 中 CR 为 压缩 比 ，PSNR. 为 峰值 信 
噪 比 . 
由 表 1 可 以 看 出 ， 在 相同 峰值 信 品 比 的 情况 下 ， 本 文 所 给 方法 的 压缩 比 要 高 于 Davis 方 法 . 


表 1: ”本文 方 法 与 Davis 编码 方法 的 结果 比较 





编码 方法 Bit Rate (bpp) CR PSNR (dB) 
本 文 方法 0.118 67.3:1 29.9 
Davis 方法 0.121 65.6:1 29.9 








图 3: 512 x 512 Lena 原始 图 像 图 4: 本 文 方法 实验 结果 CR=67.3, PSNR-29.91 


4.2. 各 类 小 波 树 统 计 

我 们 就 不 同 的 阐 值 TT， 统计 了 在 一 次 编码 中 每 一 分 辨 率 层 内 各 类 树 的 数量 ， 结 果 列 于 
表 2， 其 中 工 .为 匹配 误差 的 阔 值 ， 为 判断 小 波 系 数 是 否 有 效 的 阔 值 ，J 为 子 树 根 结 点 所 在 的 
分 解 层 . 


表 2: 各 类 树 统计 结果 ， 其 中 J 为 子 树 根 结 点 所 在 的 层 





Lena Te =50, J=5 T=40, J=4 Te=40, J=3 总 计 





Sep 167 291 1070 1528 
151 262 1031 1444 

根 重 要 树 93 202 762 1057 
87 179 770 1036 

可 自 变换 树 172 202 220 666 
192 237 299 807 

复杂 树 336 649 544 1457 


338 674 596 1529 
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由 表 2 可见 ， 在 需要 编码 的 子 树 中 根 重要 树 占 有 较 大 的 比例 ， 平 均 在 20% 左右 .这 说 明 在 
分 类 中 增加 根 重 要 树 的 分 类 是 合理 的 ， 因 为 对 根 重要 树 ， 不 需要 搜索 匹配 父 树 ， 对 其 单独 编码 
可 减少 计算 量 ， 进 而 提高 编码 速度 . 

表 3 给 出 了 每 一 层 的 搜索 次 数 (每 棵 子 树 在 父 树 库 中 搜索 一 遍 算 一 次 )， 在 采用 本 文 所 给 出 
的 方法 时 ， 搜 索 次 数 明 显 减 少 ， 如 在 低层 J = 3 时 ， 搜 索 次 数 减 少 49%. 


表 3: 搜索 次 数 统计 (RB To = 50)， 其 中 J 为 子 树 所 在 的 层 


Lena T, =50, J=5 Te=40, J=4 Te=40, J=3 总计 

本 文 方法 508 851 764 2123 
分 形 方法 601 1053 1526 3180 
减少 搜索 次 数 1596 1996 4996 3396 








4.3 匹配 映射 效果 实验 

为 了 验证 双 因 子 匹 配 映 射 的 有 效 性 ， 我 们 统计 了 分 别 采 用 双 因 子 匹 配 映射 与 单 因子 匹配 映 
射 编码 时 ， 每 一 分 辩 率 层 内 可 自 变换 树 的 数量 ， 结 果 列 于 表 4. 

由 表 4 可 见 ， 双 因子 匹配 映射 在 较 低频 子 带 可 以 找到 较 多 的 可 自 变换 树 . 因为 根 结 点 在 较 
低频 子 带 的 子 树 包含 较 多 的 小 波 系 数 ， 所 以 使 用 双 因 子 匹配 映射 具有 明显 的 优势 . 








表 4: 每 一 子 带 内 可 自 变 换 树 数 量 统计 ， 表 中 左边 Tb = 50， 右 边 T = 40 











Lena T=50 J=5 T=40 J=4 T=40 J=3 
双 因 子 匹配 映射 172 192 202 237 220 299 
单 因子 匹配 映射 120 140 143 184 123 180 





5 ”结论 


本 文 主要 结合 小 波 嵌 入 式 零 树 编 码 思 想 研 究 了 分 形 图 像 压 缩 方法 ， 提 出 了 新 的 小 波 树 分 类 
方法 与 具有 双 伸 缩 因 子 的 匹配 映射 ， 并 用 数值 实验 验证 了 所 给 方法 的 有 效 性 ， 进一步 的 工作 将 
从 理论 上 分 析 算 法 的 复杂 性 ， 这 对 于 匹配 映射 的 研究 具有 重要 意义 ， 
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Wavelet Fractal Image Compression Based on Mapping with Dual 
Scale Factors 


HUI Cun-yang 


(Department of Mathematics, Xi'an University of Arts and Science, Xi'an 710065) 


Abstract: The paper considers the fractal image compression system based on the wavelet transform. 
By the energy distribution of wavelet coefficients for images, we give a new kind of classification methods 
for wavelet trees and establish a compression mapping with two scale factors which is used for image 
coding. Using the idea of the embedded image coding using zero-trees of wavelet coefficients, a new 
wavelet fractal image compression method is proposed. Several numerical experiments are conducted 
and compared with existing methods. 
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